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ABSTRACT 

We present the first ground-based detection of sodium absorption in the trans- 
mission spectrum of an extrasolar planet. Absorption due to the atmosphere 
of the extrasolar planet HD 189733b is detected in both lines of the Na I dou- 
blet. High spectral resolution observations were taken of eleven transits with the 
High Resolution Spectrograph (HRS) on the 9.2 meter Hobby-Eberly Telescope 
(HET). The Na I absorption in the transmission spectrum due to HD 189733b is 
(—67.2 ± 20.7) X 10~^ deeper in the "narrow" spectral band that encompasses 
both lines relative to adjacent bands. The Icr error includes both random and 
systematic errors, and the detection is >3cr. This amount of relative absorp- 
tion in Na I for HD1897 33b is ~3x larger than detected for HD209458b by 



Charbonneau et al.l (120021 ) . and indicates these two hot-Jupiters may have sig- 



nificantly different atmospheric properties. 

Subject headings: atmospheric effects — line: profiles — planetary systems — 
planets and satellites: individual (HD189733b) — stars: individual (HD189733) 
— techniques: spectroscopic 



1. Introduction 

High signal-to- noise (S/N) and high spectral resolution observations of transiting ex- 
oplanets provide an opportunity to measure the properties of exoplanet atmospheres and 
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exospheres. Early theoretical models predicted that absorption, due to exoplanetary atmo- 
spheres, should be detectable in the strong resonance lines of Na I an d other alkali metals 



i n opt ical transmission spectra (ISeager &: Sasselovl I2OOOI : iBrownl I2OOII ) . ICharbonneau et al. 
( I2OO2I ) used medium resolution spectra (i? = A/AA ~ 5,540) of Na I from the Space Telescope 



Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST) to make the first 
detection of absorption in the transmission spectrum caused by the exoplanetary atmosphere 
of HD209458b during transit. The measured Na I absorption, made in the "narrow" spectral 
window (5887-5899 A) since the doublet is not fully resolved, was (—23.2 ± 5.7) x 10 ^ deeper 
than in adjacent bands during transit, and weaker than model predictions by a factor of ~3. 
Adjustments to cloud levels, metallicity, rainout of condensates, distribution of stellar flux, 
and photoion ization of sodium could account fo r the discrepancy between the models and 
observations (IBarmaru 120071 : iFortney et al.l 120031). Exosph e ric absorpti on was also detected 
for HD209458b in several UV lines bv lVidal-Madiar et al.l fl2003l . l2004h . 



Previous ground-based attempts to detect absorption in the optical transmission spectra 



due to limited observations of only a single effective transit (e.g., 


Bundv & Marcv 


2000; 


Moutou et al. 


2001; 


Winn et al. 


2004; 


Narita et al. 


2005; 


Bozorgnia et al. 


2006). In addition. 



contaminants, such as telluric lines and, in some cases, the use of an iodine cell, hindered or 
prevented measurements of several spectral features, most notably the Na I lines. 



Bouchy et al.l (120051 ) announced the identification of the short-period (P = 2.2 days) 
transiting exoplanet around the bright {V = 7.7) KOV s tar, HD189733. Subsequent work 



on this system tightly c onstrained the orbital parameters (IHebrard fc Lecavelier Pes Etangs 
20061 ; iBakos et al.ll2006l ) , determined that the axis of st ellar rotation is al igned with the plan- 
etary orbital axis using the Rossiter-McLaughlin effect (jWinn et al.ll2006l ). fo und evidence for 
photometric stellar vari ability and a stellar rotational period of ~12 da ys (IHenrv fc Winn 
20071 ; IWinn et al.l 120071 ). and identified transit signatures of starspots (jPont et al.l 120071 ). 
Infrared studies of thi s system have been successful in character izing the atmospheric tem- 
perature distribution (jPeming et al.ll2006l ; iKnutson et al.ll2007bl ). 



2. Observations and Data Reduction 



Based on the previous ground-based attempts to detect exoplanetary absorption in 
transmission spectra, our program strategy was designed to observe multiple transits and to 
specifically target the strong Na I lines. Observations of HD189733 were taken over the course 
of a year (11 Aug 20 06 to 11 Aug 2007) using the 9.2 meter Hobby-Eb erly Telescope (HET; 
Ramsey et al.l Il998l ) High Resolution Spectrograph (HRS; iTullI Il998l ). These observations 
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cover the spectral range from 5000-9000 A, with a resolving power of i? ~ 60,000. No iodine 
cell was used. A bright, rapidly rotating, B star [a Del) was observed immediately after 
HD 189733 in order to characterize the telluric absorption. We obtained spectra during 11 
in-transit visits, and 25 out-of-transit visits. Each visit is typically comprised of several 
(~3.2 during transits and ~5.1 outside of transits) 600 second exposures, and the S/N at 
Na I of each exposure is ~320 per resolution element. Observations are considered in-transit, 
if the entire exposure is obtained within first and fourth contact. 



The data were reduced using Image Reduction and Analysis Facility (IRAF; lTodylll993l ) 
and Interactive Data Language (IDL) routines to subtract the bias, flat field the images, re- 
move scattered light and cosmic ray contamination, extract the echelle orders, calibrate the 
wavelength solution, and convert to heliocentric velocities. Local interstellar medium (LISM) 
and stellar lines were removed from the telluric standard before being applied to the target 
spectra. LISM Na I absorption is detected in the spectra of a Del ( ~74 pc), our telluric stan- 
dard, consistent with LISM detections toward other nearby stars (IRedfield fc Linskvl 120021) 



and p revious Ca II and Na I observations of a Del by IVallerga et all Jl993h and lWelsh et al.l 
( 1l99ll ). Stellar lines in our telluric stand ard, a rapidly rotating hot B-star, are very weak 



and were removed using a stellar model (IKoesterke et al.l 120071 : iFitzpatrick fc Massal 12005 



Castelli &: Kuruca l2006l ). Observations of HD189733 were divided by the cleaned telluric 
standard spectra to remove the contaminating telluric features. 

To create the in-transit and out-of-transit templates, additional processing was required 
to match the blaze response and wavelength calibration of all observations. Over the short 
spectral windows (~30 A) of interest, fluctuations in the blaze response are smooth and 
slowly varying, and were easily characterized by a low order polynomial. Likewise, precise 
spectral alignment required small adjustments (~0. 1-0.2 pixels) to the wavelength solution 
based on the cross correlation of nearby stellar lines. 



3. Na I Transmission Absorption 

Figure [1] shows the observed Na I region of the stellar spectrum. The doublet is fully 
resolved and accompanied by many weak stellar lines. The bottom plot shows the transmis- 
sion spectrum, or the difference of the relative fluxes, (JTj^ — J^o^t) /J-'out, of the ultra high 
S/N in-transit template {J^in, S/N 1600) and the even higher S/N out-of-transit template 
(^out; S/N r-^ 3400). The signal predicted due to differential stellar limb- darkening is shown, 
and has been removed from the observations, although the impact is minimal. 

Excess absorption is clearly detected in the in-transit spectrum (^in) compared to the 
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out-of-tran sit spectrum (^mit) for bot h Na I lines. Using the same "narrow" spectral region 
defined by ICharbonneau et al.l (120021 ). 5887-5899 A, we measure an excess of absorption 
(—67.2 ± 7.2) X 10~^, when compared to the adj acent spectrum. This is ^ 3x larger than 



measured for HD209458b ([-23.2 ±5.7] x lO'^; ICharbonneau et all 120021 ). We measure a 
random error of 7.2 x 10~^, comparable to that achieved by the STIS/i/S'T observations. 

We also analyze a neighboring strong stellar line that is predicted to show no significant 
exoplanetary atmospheric absorption. In Figure [21 we present the same analysis near the Ca I 
line at 6122 A. No absorption is detected in this region, and over a narrow spectral range from 
6119-6125 A, we measure an insignificant amount of emission (+15.4 ± 6.8) x 10 ^, consistent 
with there being no difference between the in-transit and out-of-transit templates. 



3.1. Stellar Limb Darkening 



Differential limb-darkening, the relative difference in the limb-darkening response of the 
cores of spectral lines and the continuum, contributes to the transmission spectrum. Since we 
are dealing with a narrow spectral window, and the spectra are normalized by the adjacent 
continuum, these obser vations are not sensitive to the broadban d color dependency expected 
due to limb- darkening (IBrown et al.l 1200 ll : iBallester et al.l 120071 ). 



Esse ntially, the in-tr ansit template is a collection of spectra where a small fraction, 
(~2.5%; iPont et al.l 120071 ). of the stellar disk is blocked. This blocked component causes 
a loss of flux at a particular radial velocity, and leads to the Rossiter-McLaughlin effect 
( IWinn et al.ll2006l ). In addition, the blocked component is associated with a speciflc fi- 
angle, and therefore causes the loss of a specifically limb-darkened fraction of the stellar 
disk. The observational distribution of angles, all necessarily <0.73 sin ce the HD 189733 b 
transits at an angle slightly inclined from precisely edge-on {i = 85.7°; iPont et al.ll2007l ). 
will produce a unique differential limb- darkening signature. Indeed, with sufficient S/N, a 
transiting exoplanet could be used as a blocking obstruction, where the resulting differential 
limb-darkening spectrum would constrain, among other things, the temperature profile of 
the stellar atmosphere. In our case, we are concerned about the contribution this signal may 
have on the absorption feature in the transmission spectrum. 

We use a stellar mo del of HD 189733, ba sed on literature values. Teg = 5050 K, logg = 
4.5, and Fe/H= —0.03 ( Bouchv et al. 20051), to calc ulate the stellar fiux blocked by the 
exoplanet for each observation (jKoesterke et al.ll2007l ). We assume a 0.755 Rq star and a 
1.154 i?jup planet (jPont et al.l 120071 ). The blocked fiux is subtracted from the integrated 
disk fiux and the continuum is normalized. The colored spectrum in the bottom panel of 
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the left plot in Figures [T] and [2l shows the relative difference spectrum due to differential 
limb-darkening for our distribution of in-transit observations. The contribution of differential 
stellar limb-darkening is m uch smaller than the observed absorption, and as for HD209458 
( ICharbonneau et al.ll2002l ). at Na I, it is actually a net emission, where it contributes +1.46 x 
10~^ in the "narrow" spectral window. At Ca I the contribution is —0.10 x 10~^. In both 
cases, the differential limb-darkening contribution is well below the la random error. 



3.2. Empirical Monte Carlo 

Given the weakness of any atmospheric absorption in a transmission spectrum, it is 
critical to assess the contribution of systematic errors. Systematic errors associated with 
data reduction may arise from the continuum normalization, wavelength calibration, removal 
of the stellar and interstellar lines from the telluric standard, and the differential limb- 
darkening calculation. In addition, there are systematic contributions that are unrelated to 
data reduction procedures, including stellar variability, the changing distribution of starspots, 
and transits over active regions. These issues are particularly important since HD189733 is 
a moderately active star that is k r iown to exhibit photometric variability due to star spots 



(iHenry fc Wind 120071 : iPont et al.l 120071 ). Furthermore, because our dataset is comprised 
of many observations, spread over a year, it is important to quantify the stability of the 
individual observations that contribute to our in-transit and out-of-transit templates. 

We have constructed a diagnostic that we call an empirical Monte Carlo, or bootstrap- 
ping, analysis. It involves randomly selecting from the ~200 individual exposures a test 
sample of "in-transit" exposures, running through our entire data analysis algorithm, and 
measuring any signal in the transmission spectrum. We repeat this exercise thousands of 
times in order to establish statistical significance. Three different scenarios are explored: (1) 
an "out-out" comparison, where a subset of out-of-transit exposures, equal to the number 
of our in-transit exposures (~36), is selected and compared to the rest of the out-of-transit 
exposures; (2) an "in- in" comparison, where a subset of in-transit exposures (~8), set at the 
same ratio as the number of in-transit exposures to the number of out-of-transit exposures 
(36/165), is compared to the rest of the in-transit exposures; and (3) an "in-out" comparison, 
where an increasing number of in-transit exposures (up to ~half of the sample) is removed 
from the in-transit template and compared to the out-of-transit template. The distribution 
of the "narrow" band difference of the relative flux measurements for all realizations of the 
three scenarios are shown for Na I in the right panel of Figure [1] and for Ca I in Figure [2j 



Both the "in-in" and "out-out" distributions are centered on zero net signal, indicating 
that there is no significant difference between the different out-of-transit exposures ("out- 
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out" scenario), and no significant difference between the various in-transit exposures ("in- in" 
scenario). Since we have fewer in-transit exposures in order to do an "in-in" comparison, the 
S/N is greatly reduced, and the distribution is much wider than the "out-out" distribution. 
For Na I, the mean value (/x) of the "out-out" distribution is —3.5 x 10^^ with a standard 
deviation (cr) of 20.7 x 10~^, whereas the "in-in" distribution n = 3.9 x 10~^ and a = 
35.5 x 10^^. The results are similar for the Ca I line where for the "out-out" distribution 
/i = —0.5 X 10~^ and a = 13.7 x 10"^, and the "in-in" distribution fj, = —4.0 x 10"^ and 
a = 32.8 X 10~^. The Ca I line is located more centrally on the blaze and therefore is acquired 
with slightly higher S/N, which translates into a narrower distribution of measurements. 

As expected, the "in-out" distributions are centered on the value we obtained from the 
full dataset. The width of this distribution provides a consistency check that a handful of 
in-transit exposures are not dominating the measurement. This demonstrates, along with 
the "in-in" distribution being centered at zero, that all the in-transit exposures, over 11 
transits, consistently show excess absorption from the transiting exoplanetary atmosphere. 

The width of the "out-out" distribution measures directly the propagation of systematic 
errors through our analysis. These realizations include the contribution of both reduction 
systematics, as well as any astrophysical systematics which are not correlated with the 
exoplanet orbital period. For example, if changing distributions of starspots (i.e., different 
phases of the stellar variability cycle) could mimic an absorption feature in the transmission 
spectrum, we would expect to see a significant fraction of "out-out" realizations result in such 
measurements. Because no such realizations are able to replicate the level of the detection, 
and because the orbi tal period is completely decoupled from the stellar variability period 



(IHenry fc Wind 120071 ) such that our 11 in-transit visits are randomly distributed in stellar 
rotational phase, stellar variability cannot be responsible for the observed absorption in the 
transmission spectrum. The standard deviation of our "out-out" distribution is adopted as 
the appropriate la error for our transmission spectrum measurements since it encompasses 
both random and systematic errors. Our Na I absorption detection is (—67.2 ± 20.7) x 10~^, 
while the signal at Ca I is (+15.4 ± 13.7) x 10"^ 



4. Properties of the HD189733b Atmosphere 

The absorption feature in the transmission spectrum in Figure [T] is a direct measurement 
of the effective radius of the exoplanetary atmosphere due to the high opacity at the center 
of the strong Na I lines. The right axis of Figure [1] shows the magnitude of the absorption in 
units of exoplanetary radius. The cores of the Na I lines extend to effective radii of ~1.06 Rp\. 
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The spectra have not been corrected for the Doppler shift of the exoplanetary absorption. 
During conjunction, (e.g., mid-transit), the radial velocity of the planet is kms~^, although 
the rate of change in radial velocity is at its maximum. Over the duration of the transit, 
the radial velocity of HD189733b will change by ±19.2 kms~^. This will cause a smearing 
of the absorption profile. The observed feature appears to be significantly blueshifted from 
the stellar line center, by ~38 km s~^, and may be a combination of the planetary orbital 
motion and high speed wi nds flowing from the hot dayside, predicted by several dynamical 
atmospheric models (e.g.. lCho et al.ll2003l : ICooper fc Showmanll2005l ). 



Recent infrared observations of HD209458b by lKnutson et al.l (l2007al) show evidence for 
a temperature inversion in the upper atmosphere of that exoplanet. iBurrows et al.l (120071 ) 
presented accompanying models which require an additional absorber in the upper atmo- 
sphere of HD209458b, which could be high-altitude clouds and explain the weaker than 



expected Na I abso rption signal in the transmission spectrum (jCharbonneau et al.l 12002 



Fortney et al.ll2003l ). Our measurement of strong Na I absorption argues against the pres- 
ence of such clouds high in the atmosphere o f HD189733b, consi stent with the isothermal 
temperature profile presented for HD189733b (iTinetti et al.l 120071 ). 



5. Conclusions 



We present the first ground-based detection of absorption in the transmission spec- 
trum of a transiting extrasolar planet. The detection is based on the first high resolution 
{R ~ 60,000) transmission spectrum, with the Na I doublet fully resolved. The Na I trans- 
mission spectrum absorption for HD189733b, measured in a "narrow" spectral band that 
encompasses both lines, is (—67.2 ± 20.7) x 10~^. The estimated error includes both random 
and systematic contributions. The amo unt of absorption in Na I du e to HD189733b is ~3x 
larger than detected for HD209458b by ICharbonneau et al.l (|2002| ). and indicates that the 
two exoplanets may have significantly different atmospheric properties. The excess absorp- 
tion appears in a blue shifted component, while the red half of the Na I doublet lines has a net 
deficit of absorption during transit. Future work will include the analysis of other important 
optical lines (e.g., K I, Ha), detailed analysis of the absorption profile and comparison with 
atmospheric models, measurements of spectral indicators of stellar variability, and observa- 
tions of other transiting systems. High resolution space-based spectroscopic observations will 
provide a reliable confirmation of this detection, and access to important lines in the ultra- 
violet. Transmission spectra do not only have the potential to probe the physical properties 
of hot- Jupiter and hot-Neptune exoplanets, but lay the foun dation toward detailed physical 
diagnostics of Earth-sized extrasolar planetary atmospheres ( lEhrenreich et al.ll2006l ). 
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Fig. 1. — top left: Spectrum of HD 189733 near the Na I doublet, bottom left: The difference 
of the relative flux of the in-transit template (^in) and the out-of-transit template (^out)- 
Excess absorption in the in-transit template is clearly detected in both Na I lines. Also 
shown is the contribution of differential limb darkening (red line and filled circles), which 
has been removed from the data, but the impact is minimal. The vertical dashed lines indi- 
cate the narrow bandpass used t o calculate the relat i ve tra nsmission spectrum absorption, 
and is identical to that used by ICharbonneau et al.l (120021 ). The righthand axis gives the 
effective radius of the planet as a function of wavelength, assuming stellar and planetary 
properties from lPont et al.l (120071 ). right: Distributions of the empirical Monte Carlo analy- 
sis, which demonstrate the stability of the observations and the contribution of systematic 
errors. The full dataset transmission spectrum absorption measurement is indicated by the 
vertical dashed line. Given the random and systematic errors measured by the standard 
deviation of the "out-out" distribution, this exoplanetary absorption detection is >3a. 
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Fig. 2. — Same as Figure [T], but for the Ca I line, which hke Na I is a strong hne, but 
predicted not to show absorption in the transmission spectrum. 



